Mon. Not. R. Astron. Soc. 000,[TJ{2]() Printed 4 February 2008 (MN OTpX style file v2.2) 



The iron abundance of the Magellanic Bridge * 



oo 
o 
o 

(N 



Oh 
6 



> 



O 

OO 

O 



P.L. Dufton 1 !, R.S.I. Ryans 1 , H.M.A. Thompson 1 and R.A. Street 2 3 

1 Department of Physics and Astronomy, Queen's University of Belfast, Belfast BT7 INN, United Kingdom. 

2 Las Cumbres Observatory, 6740B Cortona Drive, Goleta, CA 93117, USA. 

3 Department of Physics, Broida Hall, University of California, Santa Barbara, CA 93106-9530, USA. 



Accepted Received ; in original form 



ABSTRACT 

High-resolution HST ultra-violet spectra for five B-type stars in the Magellanic Bridge 
and in the Large and Small Magellanic Clouds have been analysed to estimate their iron 
abundances. Those for the Clouds are lower than estimates obtained from late-type stars or 
the optical lines in B-type stars by approximately 0.5 dex. This may be due to systematic 
errors possibly arising from non-LTE effects or from errors in the atomic data as similar low 
Fe abundances having previously been reported from the analysis of the ultra-violet spectra of 
Galactic early-type stars. The iron abundance estimates for all three Bridge targets appear to be 
significantly lower than those found for the SMC and LMC by approximately -0.5 dex and -0.8 
dex respectively and these differential results should not be affected by any systematic errors 
present in the absolute abundance estimates. These differential iron abundance estimates are 
consistent with the underabundances for C, N, O, Mg and Si of approximately -1.1 dex relative 
to our Galaxy previously found in our Bridge targets. The implications of these very low metal 
abundances for the Magellanic Bridge are discussed in terms of metal deficient material being 
stripped from the SMC. 
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1 INTRODUCTION 

The Large (LMC) and Small (SMC) Magellanic Clouds are ir- 
regularly shaped galaxies, which are gravitationally bound to our 
Milky Way. In addition to these stellar systems, there are two re- 
gions of neutral hydrogen, designated the Magellanic 'Stream' and 
'Bridge'. The former traces out a large arc in the sky from between 
the C louds to beyond the South Galactic Pole (Fujimoto & Murai 
1 1984b . The latter is a b and of material which j oins t he two Clouds 
(Hind man et al. 1 19631 : Staveley-Smith et al. 1 19981 : Muller et al. 
120041) and is widely considered to be a remnant of a close tidal 
encounter of the two Cloud s abo ut ~200 Myr ago (Murai and Fuji- 
moto ll980l : Gardiner et al. 1994). This p roduced a tidal bridge and 
tail st ructure (Gardiner & Noguchi ll996l : Sawa, Fujimoto & Kumai 
1999), which is viewed as overlapping in the sky, with the tail ma- 
terial modelled to be more distant than the bridge and with a larger 
radial velocity. 
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NASA contract NAS 5-26555. These observations are associated with pro- 
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As H I emission and stellar condensation are strongly corre- 
lated, it is natural to suppose that star-formation may take place 
in these regions. Additionally, molecular clouds have been iden- 
tified in the Magellanic Bridge from CO emission (Mizuno et al. 
2006), which given the short lifetimes of such clouds were proba- 
bly formed after the tidal encounter. Attempts to find a stellar pop- 
ulation as sociat ed with the Stream have been largely unsuccess- 
ful (Irwin ll99 If) . By contrast, young stellar assoc iations have been 
identi fied in the Bridge (Irwin et al.[l985; 1990; Demers & Irwin 
Il99ll : Bica and Schmitt ll9951) and may reflect an episode of star 
formation triggered by the close encounter of the two Clouds. 

The chemical composition of both the SMC and LMC have 
been extensively studied using both stellar targets and emission 
nebulae. For example the spectra of early -type stars have bee n used , 
viz, O-type supergiants (Crowther et al. 2003; Hillier et al. 120031) . 
O-type main sequen ce tar gets (Bouret et al. 120031) . B-type super- 
giants (Trundle et al. 2004 ; Dufton et al. 120041) . B-type giants (Korn 
et al. 200d ; Lenn on et al. ll999h B-type ma in sequence stars (Korn 
et al.[2002; 2005, Hunter et al l2004l;l2007h and A-type supergiants 
(Venn 2000; Venn & Prvzbilla 120031) . Late-type stellar investiga- 
tions have included FG -type super giants (Russell & Dopita 1992; 
Spite et al. ll989l : [l99~ll; Hi ll et al. ll995l : Andrievsky et al. 120011) . 
Ceph eids (Luck et al [l998f) and K-type super giants (Barbuy et al. 
1 199 it Hill ll997l ; [T999l ; Gonzalez & Wallerstein Tl999f) . Additionally 
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th ere hav e been many H II region studies (see, for example, Kurt e t 
al. 1 19991 : Peimbert et al. |200Ct Testor l200ll ; Tsamis et al. 120031) . 
In general terms, these imply that the LMC and SMC are under- 
abundant in metals (compared with our Galaxy) by approximately 
-0.3 and -0.6 dex, with nitrogen showing a greater underabundance 
particularly for unevolved stellar targets and H II regions. 

The chemical composition of the Magellanic Bridge is more 
uncertain. Hambly et al. J 19941) obtained optical spectra for two 
early-type stars (DI 1194 and DI 1388) in order to investigate 
whether their chemical compositions reflected those of the LMC 
or SMC. However, it was only possible to make marginal identifi- 
cations of a few metal lines, and these implied large but uncertain 
metal depletions. Rolleston et al. ( 1999) analysed high and medium 
resolution spectroscopy of three other B-type bridge stars. They 
found an underabundance in the light metals (C, N, O, Mg, Si) of 
approximately -1.1 dex compared with our Galaxy. Additionally 
they undertook a differential analysis relative to the SMC B-type 
main sequence star, AV 304, and found the Bridge material to be 
relatively underabundant in metals by -0.6 dex. Unfortunately no 
iron features were observed in these optical spectra thereby pre- 
venting an abundance determination for this element. Lehner et al. 
(2001) have studied the interstellar ultraviolet absorption and H I 
emission spectra towards a young star in the Bridge. These obser- 
vations also were consistent with a low metallicity ( ~ —1. 1 dex) 
compared to that of our Galaxy. Recently Lee et al. (2005) stud- 
ied B-type targets in the SMC south east wing which adjoins the 
Bridge. These were found to have a normal SMC chemical compo- 
sition, although a re-analysis of the extant Bridge observational us- 
ing sophisticated non-LTE techniques confirmed the very low metal 
abundances found previously by Rolleston et al. 

The low metal abundance deduced for the Bridge material is 
surprising as one would expect its metallicity to reflect that of its 
progenitor material. A possible explanation is that the inter-cloud 
material could have been stripped off from either of the galaxies at 
a much earlier encounter (about ~8 Gy r ago; Kobulnicky & Skill- 
man 1 19971 : Da Costa & Hatzidimitriou 1 19981) . when less nucleo- 
synthetic processing of interstellar medium of the SMC had oc- 
curred. However, this is inconsistent wit h the numerical simula- 
tions of, f or exa mple, Murai and Fujimoto dl98d) and Gardiner and 
Noguchi d 19961) that imply that the Magellanic Bridge was formed 
relatively recently and also with the presence of early-type stars, 
implying recent star formation. 

Here we present and analyse Hubble Space Telescope (HST) 
ultra-violet spectroscopy of three Bridge B-type stars, together 
with similar targets in the SMC and LMC. The observational data 
cover the wavelength range between approximately 1900 to 2000A, 
where the spectrum is dominated by absorption lines due to iron 
group elements. Hence the principle aim of this analysis is to in- 
vestigate whether the very low light metal abundances inferred for 
the Bridge material is replicated in the heavier metals. 



2 OBSERVATIONS AND DATA REDUCTION 

HST spectroscopy of the Magellanic Bridge stars have been ob- 
tained during two campaigns that utilised the GHRS and the STIS 
spectrographs. The former obtained spectroscopy of a Bridge tar- 
get, DI 1388, in the wavelength range 1888-1930A with a spectral 
resolution of approximately 14 km s _1 whilst the latter observed 
two Bridge targets and the LMC star, LH 10-3270 between 1888 
and 1970A with a spectral resolution of approximately 15 km s - . 



These observations were supplemented with archive STIS observa- 
tions of the SMC star, AV 304, that had the same wavelength region 
and spectral resolution. The observational data are summarized in 
Table[T]with the V-magnitudes and spectral types being taken from 
Azzop ardi et al. dl975h. Gr ondin et al. (1990), Demers & Irwin 
1991 ), Demers et al. dl99 lh. Par ker et al. dl992t) . Garmany et al. 



19871) and Rolleston et al |l999). 



For the STIS observations, multiple exposures each of approx- 
imately 3000 seconds were obtained. These were combined using 
the SCOMBINE procedure in IRAF, with tests indicating that the ac- 
tual choice of parameters was not critical. Both the combined STIS 
and the GHRS spectra were input into the data reduction package 
DIPSO (Howarth et al. 1994) and cross-correlated with theoretical 
spectra (see Sect.[3]for details) to obtain the stellar radial velocity. 
These are also listed in Table Q] and are in good agreeme nt wit h 
those found previously (Hambly et al. 1994; Rolleston et al. l2003l) . 
apart f rom the LMC target, LH 10-3270. Previously Rolleston et al. 
d2002h had found a value of 197 km s , which is smaller than the 
mean LMC radial velocity of 261km s _1 (Luks & Rohlfs dl992t) 
and the estimate of approximately 310 km s _1 found here. Hence 
this star may be a single lined spectroscopic binary. The spectra 
were corrected for these radial velocity shifts prior to further anal- 
ysis. 

It was difficult to estimate signal-to-noise (S/N) ratios due the 
rich aborption line spectra. We have used synthetic and the rela- 
tively narrow lined observational spectrum of AV 304 to aid the 
identification of suitable regions and these have been fitted by low 
order polynomials to yield the S/N ratios listed in Table Q] Given 
that as it was effectively impossible to identify true continuum re- 
gions, these values should be considered as lower limits. These esti- 
mates can be compared with those from the STIS and GHRS Expo- 
sure Time calculators which predicted S/N ratios of 25 for AV 304 
and 30 for all the other targets. These predicted values are very sen- 
sitive to the adopted reddennings, which may partially explain the 
differences with those listed in TableQ] 



3 DATA ANALYSIS 

The analysis is based on grids of non-LTE model atmospheres 
that hav e been calculated using the codes TLUSTY and S YNSPEC 
(HubenydHi; Hubeny & Lanz l 19951 : Hube ny et a l. 1 19981) . Details 
of the methods can be found in Ryans et al. (2003), while the grids 
have been discussed in more detail by Dufton et al. ( 2004). 

Briefly four grids have been generated with base metallicities 
corresponding to our galaxy ([§] = 7.5 dex) and with metallicities 
reduced by 0.3, 0.6 and 1.1 dex. These lower metallicities were 
chosen so as to be representative of the LMC, SMC and Bridge 
material. For each base metallicity, models have been calculated 
covering a range of effective temperature from 12000 to 35 000 K 
in steps of 1000K or less and logarithmic gravities (in cm s -2 ) from 
4.5 dex down to close to the Eddington limit (which will depend 
on the effective temperature) in steps of 0.25 dex. For each pair 
of effective temperatures and surface gravities, five models with 
microturbulences (£ t ) of 0, 5, 10, 20 and 30 km s~ were generated. 

Atmospheric parameters (summarized in Table O for all our 
targets have been previously estimated from optical spec tra b 
Rolleston et al. dl999l) for the Bri dge st ars, Rolleston et al. d200: 
for LH 10-3270 and Hunter et al. d2004l) for AV 304. The first two 
analyses used an LTE approach based on the model atmospheres 
of Kurucz d 1991b . whilst the last was based on the non-LTE grid 
utilised here. However atmospheric parameters for AV304 esti- 
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Table 1. Details of the HST observations. 



Star Region RA DEC Sp. Ty. V E(B-V) Inst. Proposal Wavelength Exp. Time v LSR S/N 

2000.0 mag mag A seconds km s — 1 ratio 



AV 304 


SMC 


1 02 21.4 


-73 39 15 


B0.5V 


14.98 


0.03 


STIS 


7759 


18f 


S8-1978 


5940 


145 


30 


DGIK975 


Bridge 


4 19 58.6 


-73 52 26 


B2 


15.05 


0.12 


STIS 


9130 


18* 


58-1978 


21075 


190 


29 


DI 1239 


Bridge 


2 30 40.8 


-74 04 47 


Bl 


15.24 


0.08 


STIS 


9130 


18* 


58-1978 


21075 


159 


37 


DI1388 


Bridge 


2 57 11.9 


-72 52 55 


B0 


14.39 


0.08 


GHRS 


6073 


185 


58-1929 


1410 


140 


18 


LH 10-3270 


LMC 


5 57 21.2 


-66 25 01 


Bl V 


14.90 


0.10 


STIS 


9130 


18* 


58-1978 


20480 


313 


42 



Table 2. Atmospheric parameters and projected rotational velocities 
adopted for our targets. 



Star 


T e ft 


logg 




v sin (' 




K 


cm s — 2 


km s — 1 


km s _1 


AV304 


27 500 


3.90 


3 


11 


DGIK975 


20000 


3.60 


5 


120 


DI 1239 


24 000 


3.80 


5 


90 


DI1388 


32 000 


4.00 


5 


180 


LH 10-3270 


29 500 


4.10 


6 


30 



mated by Rolleston et al. d2003h using Kurucz's LTE models are 
in excellent agreement with those of Hunter et al. and hence the 
adoption of LTE estimates should be acceptable; we return to this 
point when considering possible sources of errors in our iron abun- 
dance estimates. 



Using a grid model with an appropriate metallicity and atmo- 
spheric parameters closest to the values listed in Table [2] we have 
synthesised the spectral region observed by HST for each star. As 
the spacing between grids points is smaller than the typical uncer- 
tainties in estimating atmospheric parameters, the adoption of the 
nearest grid point is unlikely to be a significant source of error. For 
elements (carbon, nitrogen, oxygen, magnesium and silicon) ex- 
plicity included in the TLUSTY calculation, these calculations were 
in a non-LTE approximation. Iron was also included in the model 
atmosphere calculations but due to the complexity of its Grotrian 
diagram, individual levels were combined together into superlevels. 
Hence for the spectral synthesis, iron and all other elements were 
treated in LTE. Test calculation showed that typically 80% of the 
absorption in this spectral region arose from iron with over 95% 
being due to iron group elements - hence these are effectively LTE 
calculations. To fit an observed spectrum, calculations were under- 
taken for the set of element abundances adopted for a particular 
grid (see Dufton et al. 1 20041 for details) and with all these abun- 
dances scaled. The actual element abundances for any given star 
may not vary in step but as iron dominates the absorption, any dif- 
ferential variations in element abundances should not significantly 
affect these theoretical spectra. 

This approach has the consequence of introducing an incon- 
sistency between the abundances adopted in the model atmosphere 
and spectrum synthesis calculations. For example, for the star, 
AV 304, a model would be adopted from the grid with a metallicity 
appropriate to the SMC viz. -0.6 dex that of the Galaxy. Spectra 
would then be synthesised for different metallicities which might 
range, for example, from the Galactic composition to -1.2 dex of 
the Galactic composition. We have investigated the consequences 
of this simplification as folows. Firstly we calculated the spec- 



trum for an SMC model of AV 304 with atmospheric parameters, 
T e ff =27 500, log g=3.9, £t=5 km s _1 and with a metallicity -0.6 dex 
that of the Galaxy. Then we used the Galactic model with the same 
atmospheric parameters but scaled the metallicity in the SYNPLOT 
spectrum synthesis calculation in order to mimic that of the SMC. 
The resulting spectra were very similar with the summed equivalent 
widths differing by 3%. Hence we do not expect that this simpli- 
fication (th at ha s been widely used elsewhere - see, for example 
Kilian et al. ll994l Korn et al. l2000l McErlean et al l 19991) . will be a 
source of significant uncertainty. However it is possible that this ap- 
proach could lead to uncertainties that depended on the wavelength 
region observed as the blanketing will depend on the iron abun- 
dance adopted. In turn this could affect the ionization structure and 
thereby the behaviour of absorption lines formed in different parts 
of the stellar atmosphere. 

As discussed in Sect.[2] it is difficult to distinguish continuum 
regions due to the richness of the ultra-violet absorption line spec- 
trum. Hence we have not attempted any continuum normalisation 
but have instead worked directly with the flux calibrated spectra. 
For the observed spectra, we must first correct for the effect of in- 
terstellar extinction. We adopted values of the colour excess (sum- 
marized in Table[]} for the Bridge targets from Demers & Battinelli 
dl998l) . for LH 10-3270 fro m Park er et al. dl992h and for AV304 
from Azzopardi & Vigneau There is evidence that the ex- 

tinction law for the Magellanic Clouds may differ from that for the 
Galaxy (see, for example, Seaton ll979l . Howarth ll983l) . However 
given the low reddenings of our targets the use of different laws 
was found to have only a small effect on the absolute fluxes de- 
duced. More importantly, because of the limited coverage of the 
observed spectra, the relative flux distributions, which are impor- 
tant for our methodology as discussed below, remain effectively 
unchanged when corrected for extinction. Hence uncertainties in 
these corrections are unlikely to seriously affect our analysis. We 
then rebinned our unreddened spectra to a wavelength interval of 
0.1 A and to a wavelength range of 1890-1970A for the STIS spec- 
tra and 1890- 1928 A for the GHRS spectra. 

To compare our observed and theoretical spectra, we adopted 
the following methodology. The theoretical spectra were first con- 
volved with a Gaussian broadening function to allow for intrumen- 
tal broadening and with a rotational broadening function with the 
projected rotational velocities, v sin being summarized in TablefJ] 
These were taken from Hambly et al. 1 1994) for DI 1388, Rolleston 
et al. d 19991) for DI 1239 and DGI K975 , Rolleston et al. d2002ri for 
LH 10-3270 and Rolleston et al. d2003l) for AV 304. Recent stud- 
ies (s ee, for example, Dufton et al. l2006i . Simon-Diaz and Herrero 
|2007|) have shown that macroturbulence can also be important in B- 
type supergiants. The estimated gravities of our targets imply that 
they are either near to the main sequence or giants and hence that 
macroturbulence will not be significant. Additionally we use these 
projected rotationally velocity estimates to allow for the overall 
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broadening of the line spectra and to first order its physical origin 
is not important. The theoretical spectra were then rebinned to the 
same wavelength region and interval as those of the observed spec- 
tra and scaled so that they had the same mean flux as the observed 
spectra, thereby making the spectra comparable without any need 
for rectification or normalisation. The sum of the squares of the 
differences between observation and theory (S) was then computed 
and the metallicity varied until this was minimised. These metallic- 
ity estimates (designated as A[tt] because our ultra-violet spectra 
are dominated by iron group absoption lines) are summarized in 
Tabled 

This approach is illustrated in Figs. Q] and [2] for the LMC tar- 
get, LH 10-3270 and the Bridge target, DGIK975 chosen to cover 
a range of metallicities and projected rotational velocities. In the 
former the observed and best fitting theoretical spectra are shown, 
whilst in the latter, the value of S is plotted against metallicity and 
minima are found at metallicities (relative to Galactic) of -0.9 and 
-1.7 dex respectively. In general the theoretical spectra reproduce 
that observed but agreement for specific features is variable. The 
absoption line spectrum for LH 10-3270 is relatively well observed 
due to both its low projected rotational velocity and the relatively 
high metallicity of the LMC. The theoretical spectrum is in ex- 
cellent agreement with the majority of absorption features well 
matched. This in turn provides indirect support that our uniform 
scaling of all the individual element abundances is acceptable. 

For DGIK 975, the situation is less satisfactory with the match 
to individual features being far less convincing, which may be man- 
ifested in the broad minimum found in the goodness-of-fit shown 
in Fig. [2] Given the better agreement found for LH 10-3270, we 
ascribe these problems principally to the observed absorption lines 
being intrinsically weak and broadened by the significant stellar 
rotation. However part of these discrepancies may also be due to a 
combination of inadequacies in the atomic data (and in particular 
the oscillator strengths) and the neglect of non-LTE effects for the 
species that dominate the absorption in this spectral region. 

The spectra line data used in SYNSPEC is mainly tak en from 
the linelists of Kurucz (see, for example, 1 19881 : Il99lh . These 
have been updated using more recent results taken from the 
NIST databa se as discussed in for example Lanz and Hubeny 
d2003l : |2007|) . with further information being available from the 
TLUSTY/SYNSPEC webpage. All the energy levels are experimen- 
tal and as such the corresponding wavelengths of the transistions 
should be reliable. This was confirmed by comparing the our wave- 
lengths for a sample of Fe III featur es wit h those determined ex- 
perimentally (see, Sugar and Corliss. 1 19851 and references; Ekberg, 
Il993l therein) where discrepancies were typically less than 0.01 A. 
This would correspond to a velocity uncertainty of approximately 
1 km s _1 , which is far smaller than the intrinsic width of the line 
due to the rotational broadening. Hence we do not believe that un- 
certainties in the wavelengths will be a significant source of error. 

In the case of the oscillator strengths, a number of 
indepe ndent calculations exist apart from those of Kuruc z 
dl99ll) , i ncludi ng Toner and Hibbert J2005h. Ekberg Jl993h . 
Fawcett (1989), Nahar and Pradhan Jl996h . Raassen and 
Uylings have also provided extensive listings via a web-site 
(http://www.science.uva.nl/pub/orth/iron/FeIIIl as part of the FER- 
RUM consortium (Johansson et al.. l2000l) . As discussed by Toner 
and Hibbert, discrepancies are typically 40%, which for a single 
feature would lead to an absolute abundance error of appro ximately 
0.2 dex. However as discussed by Thompson et al. j2007h . for the 
stronger lines the discrepancies are smaller, being typically 10% or 
less. Additionally as we fit a large number of features simultane- 



ously, the corresponding abundance error would be smaller if these 
discrepancies were random. The errors in differential abundances 
should also be smaller as errors due to uncertainties in the oscilla- 
tor strengths should to first order cancel. As even for the best mod- 
ern calculations, there remains differences between the theoretical 
oscillator strength estimates and as these new results cover only a 
subset of the lines incorporated in our linelists, we have decided 
not to attempt to update our linelist. Given the above discussion the 
abundance estimates (listed in Table [3]( must be treated with cau- 
tion. However variations in these estimates should be more reliable 
as they are based on a consistent and objective approach. 

The above methodology implicitly assumes that the flux cali- 
bration of the observed spectrum has been reliably performed. For 
example there could remain a low frequency ripple and to investi- 
gate the consequence of this we have artificially added a sine wave 
component with an amplitude of 1 % and a period equal to the wave- 
length range of our data for the target, LH 3270. We then repeated 
the analysis and found the same abundance estimate (after round- 
ing to the nearest 0.1 dex). Clearly larger errors would map onto 
larger uncertainities in the abundances but we would expect that 
any systematic errors in the calibration would similarly affect all 
our abundance estimates and hence relative values should be more 
reliable. 

One significant discrepancy in the comparison shown in Fig.Q] 
is that the theoretical spectra predicts a relatively strong feature at 
approximately 1906. 8 A. This is not apparant in the observed spec- 
tra for LH 3270 and also may be absent from that for DGIK 975 
although in this case the comparison is complicated by the greater 
amount of rotational broadening and the intrinsic weakness of the 
spectrum. Investigation of the theoretical spectra indicates that this 
feature arises principally from an Fe III line, which has been in- 
corporated with and oscillator strength 0.78. For this transition, 
Raassen and Uylings derive a smaller value of 0.41, which may 
at least partially account for the discrepancy. As a test, we have re- 
peated the fitting procedure for LH10 3270 excluding this region. 
The goodness-of-fit was indistinquishable from that shown in Fig. 
[2] and lead to the same estimate of the iron abundance as listed in 
Tabled 

For DI 1388, the wavelength coverage of its GHRS spectrum 
is approximately half of that of the STIS observations and this leads 
to its metallicity estimate not being strictly comparable with those 
for the other targets. Additionally the S/N ratio of the GHRS spec- 
trum is relatively low, whilst DI 1388 has the highest vsin/ esti- 
mate. Hence the metallicity estimate for this star is particularly 
uncertain, although we note that from visual inspection the agree- 
ment between observation and theory appears poor for values of the 
metallicity, A[-^] > -1.3 dex. 

Inspection of Table [2] indicates that the two SMC and LMC 
targets have smaller projected rotational velocities and hence their 
spectral lines will be narrower. In turn this could affect the compar- 
ison between iron abundances deduced in these targets and those 
in the Bridge. We have investigated this as follows. For these two 
targets we have rotationally broadened their observed spectra to a 
value of vsin/= 120 km s , which is typical of that observed in 
our Bridge targets. We have then added additional Poisson noise to 
obtain a S/N ratio of 30 again typical of the Bridge targets. We have 
then repeated the analysis of these modified spectra to obtain esti- 
mates of the iron abundances. Although there are small changes in 
the position of the minimum obtained in the goodness-of-fit, after 
rounding the abundance estiamtes to the nearest 0.1 dex, we ob- 
tain the values listed in Table [3] Hence we do not believe that the 
differences in broadening in the observed spectra should affect es- 
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Figure 1. Comparison of observed and theoretical spectra for two of our targets. The upper panel shows the observed (solid line) and the theoretical (dotted 
line) spectra for LH103270 and a metallicity of -0.9 dex (corresponding to the best fit between observation and theory) relative to that of the Galaxy. To aid 
comparison the theoretical spectrum has been shifted downwards on the y-axis. The lower panel shows the corresonding spectra for DGIK975, a metallicity 
of -1.7 dex now being adopted for the theoretical calculations 



timates of the abundance differences between the Bridge and the 
Cloud targets. 



We have estimated the uncertainties in our Fe abundance es- 
timates arising from uncertainties in the adopted atmospheric pa- 
rameters. We assumed errors of ± 1 000 K in T e s , ±0.2 dex in log g 
and ±3 km s _1 in £ t , which is consistent with the error estimates 
discussed in the analyses summarized in Table [2] We have used 
a relatively simple procedure by first calculating the total equiv- 
alent width for the ultra-violet spectra region that we have anal- 
ysed for our adopted stellar atmospheric parameters and our esti- 
mated iron abundances (see Tables|2]and[3]respectively). We then 
increased in turn each of the atmospheric parameters by its uncer- 
tainty and changed the iron abundance until we recovered the same 
total equivalent width. Tests comparing this methodology with that 
of changing each atmospheric parameter and then recalculating the 
goodness-of-fit between theory and observation for different Fe 
abundances showed that the two methods yielded similar error esti- 
mates. Our error estimates are summarized in Table[4] with values 
for a similar decrease in the atmospheric parameters being of simi- 
lar magnitude but with an opposite sign. In general, the errors aris- 
ing from uncertainties in the surface gravity appear not to be sig- 
nificant. Those for the microturbulence depend on the strength of 
the Fe absorption lines with those for the effective temperature de- 
pending on how far the effective temperature is from that (approxi- 
mately 18 000 K) corresponding to the maximum in the strength of 
the Fe spectrum. Assuming that these errors are uncorrelated, they 
can be added in quadrature and these error estimates are shown the 



final column of the Table|4] Adopting different errors estimates for 
the atmospheric parameters or allowing for the correlation between 
the effective temperature and gravity estimates would lead to dif- 
ferent estimates but we believe that a typical uncertainty of ±0.2 
dex implied by Table|4]to be reasonable. 

We have also attempted to use the value of the \ 2 associated 
with our comparison of the theoretical and observed spectra to es- 
timate the uncertainty due for example to observational errors. As 
discussed by Avni dl976l) . a change of 1.0 and 2.7 in the \ 2 pa- 
rameter should correspond to the 68 and 90% limits in our uncer- 
tainty. However utilisation of this procedure led to unrealistic small 
error estimates of less than 0.1 dex for both cases. This may re- 
flect that some of our data points (e.g. those near the continuum) 
do not represent signal but tests attempting to exclude such regions 
again led to low error estimates. Additionally the small implied er- 
rors may be due to stochastic noise not being the dominant source 
of uncertainty. For example other systematic errors in the observa- 
tional data or from the physical assumptions (including the adopted 
atomic data) may be significant. 

Moehler et al. d 19981) in their analysis of the equivalent widths 
of individual iron features in the ultra-violet spectrum of PAGB 
stars estimated errors due to observational uncertainties of 0.05 to 
0.1 dex. From inspection by eye of the quality of our fits for dif- 
ferent Fe abundances an error of approximately ±0.2 dex would 
appear reasonable which is rather larger than the errors quoted by 
Moehler et al. dl998l) . However adopting a typical random error of 
±0.2 dex and combining in quadrature with those due to uncer- 
tainties in the atmospheric parameters leads to an estimated error 
of typically ±0.3 dex. We note that this estimate does not include 
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Table 3. Estimated iron abundances relative to that of our Galaxy (i.e. A[^] = dex would correspond to a solar metallicity of 7.5 dex - Grevesse and 
Sauval 1998) for our targets deduced from the HST ultra-violet spectra. Mean light metal abundance estimates previously deduced for these targets (designated 
A[-tf] - see text for details) are also sum marized . Other estimates of the iron ab unda nces for the LMC and SMC relative to ou r Gala xy are tab ulated . These are 
from B-type supergiants by Trundle et al. 1 2007), A-type supergiants by Venn ( 2000), F-type supergiants by Andrievsky et al. 1200 II) and Hill jl997l) . The final 
colum n lists abundance studies of late-type stars, viz. Cepheids in the LMC (Luck et al.[l998) and K-type supergiants in the SMC (Gonzalez and Wallerstein 
[1991 . 
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dex 
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dex 
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dex 
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dex 
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dex 
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dex 
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-0.3 
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-0.6 
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-1.1 










DI 1388 


Bridge 


-1.9 


<-0.9 
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Table 4. Uncertainties in our estimated iron abundances due to uncertainties 
in the atmospheric parameters. For each star, the effective temperature, log- 
arithmic gravity and microturbulence have been increased by 1 000 K, 0.2 
dex and 3 km s — 1 respectively and the change of iron abundance required to 
obtain the best agreement between theory and observation was estimated; 
see text for further details. Also shown is the estimated errors assuming that 
these errors can be added in quadrature. 
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AT eff 
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A logg 
+0.2 dex 
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Quadrature 


LH 10-3270 
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-0.09 
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+0.17 
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-0.16 


±0.24 


DGIK975 


-0.06 


+0.03 


-0.12 


±0.14 


DI1239 


+0.09 


+0.03 


-0.12 


±0.15 


DI1388 


+0.18 


-0.06 


-0.02 


±0.19 



Metallicity (dex) 



Figure 2. Estimation of the metallicity for two of our targets. 
The upper panel shows the goodness-of-fit, S (with units of 10 — 28 
erg 2 cm _4 s — 2 A -2 ) between the observed and theoretical spectra for 
LH 10 3270 plotted against the metallicity adopted for the theoretical cal- 
culations. The lower panel shows the analagous results for DGIK 975 



systematic errors from, for example, the assumption of local ther- 
modynamic equilibrium in the radiative transfer calculetion and we 
return to this point in Sect. [4] 



4 DISCUSSION 

In Table[3] we also summarize for our targets, the mean differential 
abundance (designated A[t|]) for all reliably observed elements 



from carbon to sulphur with the exception of nitrogen. This is be- 
cause analysis of H I I region (see, for example, Kurt e t al. 1 19991 : 
Peimbert et al. l200d : Testor l200ll ; Tsa mis et al. 120031) and stel- 
lar spectra (see, for example, Korn et al. 120021 : 1200517 indicate that 
nitrogen has anomalously low differential abundances in the Mag- 
ellanic system. The sources of these light element abundanc es are : 
LH 10-3270 (Rolle ston et al. l2002h : AV 304 (Rolleston et al.O 
Hunter et al. l2004l) ; DKIK975 and DI 1239 (Rolleston et al. 



2003; 



1999) 



and DI 1388 (Hambly et al. 1 19941) . In this Table, we also summa- 
rize iron abundances estimated for some representative analyses of 
other stellar targets: SMC and LMC B-type super giants (Trundle 
et al.|2007); SMC A-typ e supe rgian ts (Ve nn l200rj|) ; F-type super- 
giant s (An drievsky et al. 1200 it Hill ll997t) : LMC Cepheids (Luck 
et al.[T 998) and SMC K-type supergiants (Gonzalez & Wallerstein 

A comparison of our iron abundances deduced from ultra- 
violet spectra with those obtained for other stellar targets shows 
that the former are lower by 0.5-0.6 dex for both the SMC and 
LMC. This would appear to be the case even when we compare 
the analyses of isolated optical Fe III lines of B-type stars with 
our ultra-violet spectra (which are again dominated by Fe III lines), 
although unfortunately there are no stars in common between the 
analyses. Studies of the ultra-violet spectra of other hot stars have 
also found relatively low iron abundances, as discussed by Thomp- 
son et al. d2007l) . For example, Grigsby et al. Jl998l) analysed the 
iron lines between 1800 and 2500 A in the spectrum of the B3 V 
star 1 Her. They found an underabundance of -0.47 dex relative to 
the sun assuming a zero microturbulence. As pointed out by Grib- 
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sby et al., an increase in the microturbulence would decrease their 
estimates and hence the actual underabundance estimated from the 
ultra-violet spectrum could be larger. From their results, we can 
deduce an abundance estimate from the Fe III lines that lie in the 
wavelength region (1888 to 1978 A) that we have considered. This 
again implies that iron is underabundant by 0.4 dex or more (de- 
pending on the microturbulence adopted) co mpare d with the so- 
lar value Meyer 1 19851 : Anders and Grevesse 1989; Grevesse and 
Sauval lT998h . By contrast, a nalyses of the optical Fe ill sp ectrum 
in i Her (Ko deira and Scholz ll97(ll : Pete rs and Alter ! 1970l : Peters 
and Polidan 1985; Pintado and Adelman [l985h have led to values 
that are close to solar. 

Hot post-asymptotic-giant branch (PAGB) stars have also been 
found with relatively low iron abundances as discussed by Moehler 
d200lh . Napiwotzki et al. i 1994) analysed optical and ultra-violet 
spectra of the field star BD+33°2642 and found underabundances 
(relative to solar) of approximately 1.0 dex for C, N, O, Mg, Si but 
a larger underabundance of 2.0 dex for Fe. The latter was found 
from analysis of Fe III lines in the same spectral region as consid- 
ered here. Napiwotzki et al. showed that such an abundance pat- 
tern could be understood by a seperation of gas and dust as pro- 
posed for cooler PAGB stars (see, for example, Mathis and Lamers 
l2002h . Moehler et al. Jl998h analysed ultra-violet and optical spec- 
tra of two PAGB stars, ROA 5701 and Barnard 29, in globular clus- 
ters with the main aim of determining their iron abundance. From 
the ultra-violet spectra, they determined iron abundances that for 
Barnard 29 was 0.5 dex lower than that found for the cooler giant 
stars in M 13 and for ROA 5701 was lower than the range of abun- 
dances found in uo Cen. Moehler et al. again interpreted these ob- 
servations as due to a gas/dust separation in the atmospheres of the 
PAGB stars. Recently Thompson et al. ( 12007!) have reanalysed the 
spectroscopy for these stars. From the ultraviolet spectra they find 
iron abundances using the spectra synthesis technique used here 
that are similar to those estimated by Moehler et al. However they 
also tentatively identified Fe III lines in their optical spectra. These 
implied higher iron abundances that agree well with those found 
from the red giants in M 13 and with the range of iron abundances 
found from the red giants in lu Cen (Origlia et al. l2003h . 

In summary analyses of the ultraviolet Fe III spectrum of B- 
type stars around 1900 A have consistently yielded relatively low 
iron abundances. For the highly evolved objects, this may reflect 
real physical effects but in others there appears to be differences be- 
tween abundance estimates from the optical and ultra-violet spec- 
tra and for the younger massive stars inconsistencies with the iron 
abundance that would be expected. Indeed our results for the SMC 
and LMC targets and those discussed above for i Her would imply 
that a systematic increase in the iron abundances deduced from the 
ultra-violet spectra of approximately 0.5 dex would appear appro- 
priate. The cause of this systematic offset could be in the theoretical 
modelling and in particular in the use of an LTE approximation in 
the spectral synthesis or in the adopted oscillator strengths. Addi- 
tionally if the ionic linelists are incomplete, the presence of a large 
number of weak unresolved lines would act as a pseudo-continuous 
opacity. In turn this could lead to the strengths of lines or blends 
in our theoretical spectra being too strong, thereby leading to our 
abundance estimates being too small. 

Given the uncertainties, in the absolute iron abundances de- 
duced from the HST spectra, we have instead chosen to take a dif- 
ferential approach. For example, we find a difference between in 
the iron abundance of our SMC and LMC target of -0.3 dex which 
is in excellent agreement with those found from the optical spectra 
of BAF-type and late-type stars in the Magellanic Clouds (see Ta- 



ble[3](. More importantly all three Magellanic Bridge targets appear 
to have very low iron abundances and in particularly values that are 
lower than the SMC iron abundance by between -0.4 to -0.7 dex. 
These differences are larger than the error estimates discussed in 
Sect. [3] and although it is possible that these have been underesti- 
mated, the similarity in the iron underabundance found for all three 
Bridge targets is striking. Additionally as we have adopted a differ- 
ential approach, any systematic errors should have been minimised. 

Such an iron underabundance is consistent with those for the 
lighter metals (listed in Table [3} and with the analysis discussed 
in Rolleston et al. dl999h who estimated a mean light metal under- 
abundance for the Bridge material of - 0.60±0 . 18 dex with respect 
to the SMC. Additionally Lehner et al. ( 1200 If) , from an analysis of 
the interstellar absorption lines seen in the spectra of a background 
object, also found evidence for low metal abundances in the Mag- 
ellanic Bridge gas. Therefore in the discussion that follows, we will 
assume that the material in the Magellanic Bridge is underabundant 
in metals by factors of approximately -0.5 and -0.8 dex compared 
with the SMC and LMC respectively. 

As discussed in the introduction, both o bservational (e.g. Ir- 
win et al l 19851 : fl99Cl ; Demers & Irwin ll99ll) a nd theoretical (e.g. 
Gardiner and Noguchi 1 19961 : Sawa et al. ll999h studies imply that 
the Bridge was formed during a close encounter between the Mag- 
ellanic Clouds approximately 200 million years ago. This scenario 
has been supported by more recent investigations, with for example 
Nishiyama et al. (2007) tentatively identifying pre-main sequence 
Herbig Ae/Be stars in the Magellanic Bridge. Additionally the re- 
cent photometric investigation of Harris d2007l) failed to identify an 
older stellar population implying that the material stripped from the 
Magellanic Clouds was effectively a pure gas. 

Of particular relevance to the metal abundances found in the 
Magellanic Bridge are the recent numerical simulations of Bekki 
and Chiba (2005: 120071) . These follow the evolution and star for- 
mation history of the Magellanic System over the last 800 million 
years using chemodynamical simulations with both Clouds being 
modeled as self-gravitating barred and disk systems. A mean metal 
abundance for the Bridge (relative to our Galaxy) of -0.8 dex is 
predicted, compared with -0.6 dex adopted for the cental regions of 
the SMC. This underabundance arises from the Bridge material pri- 
marily originating from the outer regions of the SMC with the sim- 
ulations incorporating a negative metallicity gradient. As pointed 
out by Bekki and Chiba <2007D . this is qualitatively consistent with 
the Br idge light element underabundances found by Rolleston et al. 
dl999h and summarized in Table [3] However as for our iron abun- 
dance estimates of the three bridge stars, there remains a significant 
quantitative discrepancy of approximately 0.3 dex. 

Bekki and Chiba d2007l) predicted a range of metallicity in 
Bridge stars with a mean value of -0.8 dex but with approximately 
three percent of stars having values ^ —0.9 dex. Our three tar- 
gets are spatially well separated across the Magellanic Bridge (see 
Figure 3 of Rolleston et al. 1 19991) and hence it would be highly 
improbable that they could all be from the lower metallicity part of 
this distribution. We therefore conclude that a plausible explanation 
for the very low light metal and iron abundances found in the Mag- 
ellanic Bridge is that the corresponding material was stripped from 
the outer regions of the SMC approximately 200 million years ago. 
However observational abundance estimates and theoretical predic- 
tions still differ by 0.3 dex. 
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5 CONCLUSIONS 

We have presented HST ultra-violet spectroscopy of five targets in 
the Magellanic Clouds and Bridge. These data have been analysed 
using grids of non-LTE models atmospheres with the spectral syn- 
thesis being predominantly undertaken in an LTE approximation 
to deduce metallicities (which will be dominated by the iron abun- 
dances). The principal results can be summarized as follows: 

(i) Iron abundances estimated from the ultra-violet spectra of 
SMC and LMC targets are typically -0.5 dex lower than those de- 
rived by other methods. This is consistent with results found from 
previous studies of Galactic objects. 

(ii) The iron abundance estimates for the three bridge targets are 
-0.5 and -0.8 dex lower than those for the SMC and LMC respec- 
tively. This is consistent with the underabundance in light metals of 
-1.1 dex found for the Bridge material with respect to our Galaxy. 

(iii) These metal abundances are in qualitative agreement with 
numerical simulations which imply that most of the Bridge mate- 
rial has been stripped from the outer layers of the SMC. However 
there remain significant quantitative differences between the ob- 
served estimates and those predicted by the simulations. 
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